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Abstract 

The polarization and the susceptibility of a ferroelectric multilayer with a non- 
ferroelectric slab are investigated within the framework of transverse Ising model 
with a four-spin interaction term. The effect of the thickness and the position of the 
non- ferroelectric slab are investigated in this paper. We find that the increase of the 
thickness of the non-ferroelectric will decrease the polarization and the susceptibility 
of the film. If the position of the non-ferroelcetric slab shifts from the center of the 
film to the surface, the number of the peaks of the susceptibility will change. And 
a step-like polarization curve is found. 
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I. INTRODUCTION 



During the past decade, there has been considerable interest in study of the ferroelectric 
heterostructure, such as ferroelectric film, multilayer, and superlattice. Heterostructure 
ferroelectric materials have important applications in ferroelectric memory, sensor, mi- 
croelectromechanical system, and optical devices. Shen et al.[l] studied a ferroelectric 
sandwich structure using the Ginzburg-Landau phenomenological theory, they investi- 
gated the effects of the long-range interaction and the thickness of the middle layer on 
the polarization of the ferroelectric sandwich structure. Wang[2,3] and Kaneyoshi[4] ap- 
plied the transverse Ising model(TIM) to investigate the phase transition properties of 
the ferroelectric film, ferroelectric/ferroelectric superlattice, and ferroelectric bilayer. The 
ferroelectric/paraelectric superlattice was studied in experiments [5-9] and in theories[10- 
11]. The enhancement of remnant polarization is found in the epitaxial BaTiO^j ' SrTiO^ 
superlattice with asymmetric structure [9], which is about three times of the BaTiO^ 
single-phase film formed under the same condition. And they proposed that the ferroelec- 
tric/paraelectric strained superlattice or multilayer was one of the effective approaches 
to enhance the remnant polarization of ferroelectric heterostructure. However, No study 
on the ferroelectric film with a non-ferroelectric slab is found in the literature. The non- 
ferroelectric slab inside the ferroelectric multilayer may be impurity or polymer, etc., 
which is inserted to the ferroelectric film to enhance the polarization and dielectric prop- 
erties, or to improve the mechanical properties of the ferroelectric multilayer. Here we 
define the inserted slab as a non-ferroelectric slab uniformly. If a non-ferroelectric mate- 
rial is introduced into the ferroelectric film, what kind of interesting phenomenon, if any, 
will happen? What is its effect on the surface and the susceptibility of the film? This 
is the motivation of our paper. As is well known, the long-range interactions dominate 
in the ferroelectric materials [12- 14]. We use the transverse Ising model with long-range 
interaction to investigate the ferroelectric film with a non-ferroelectric slab. Since most fer- 
roelectric heterostructure, studied in the experiments, are the ferroelectric material with 
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first-order phase transition, four-spin interaction [15, 18, 19] is included in the transverse 
Ising model (TIM). The TIM is simple, however it is convenient. It has been successfully 
applied in the studying of the surface and size effects in finite ferroelectric system. We 
use the TIM to investigate the effects of the thickness of the ferroelectric slab and that 
of the non-ferroelectric slab on the average polarization, the surface polarization, and the 
susceptibility of the ferroelectric multilayer. We obtain that the thicker non-ferroelectric 
slab will lower the polarization, phase transition temperature, and the susceptibility of 
the film. The position of the non-ferroelectric slab will have an important role on the 
polarization and the susceptibility of the film. A step-like polarization is found when 
the non-ferroelectric slab sites on an asymmetric position. We give a detailed discussion 
about the origin of the step-like polarization in Sec. III. 



II. Model and Formulation 

We consider a sandwich structure composed of two-sided ferroelectric slabs and a 
middle non-ferroelectric slab with different number of layers. Each layer is defined on the 
x-y plane. And pseudo spins site on a square lattice(see Fig. 1). For ferroelectric with a 
first-order phase transition, a term of four-spin interaction[18,19] should be included in the 
transverse Ising model. The system can be described by the following Ising Hamiltonian 
in the presence of a transverse field, 

h = - E Jvs;s; - £ J t3kl s?s*s z k sf - £ cusf - ^e £ s?, (1) 

<ij> ijkl i i 

where fli is the transverse field. Sf and Sf are components of spin-1/2 operator at site i. 
ji is the dipole moment on site i, and E is the applied electric field. Jy and are the 
two-spin and four-spin interaction coupling constant respectively. Since the properties 
of the surface and the interface are different from those of the bulk, we distinguish the 
interactions of the pseudo spins and the transverse field on the surface and the interface 
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from those inside the ferroelectric slabs. 
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for i,j on the surface layer, 

for i,j in slab A, 

for i,j on the interface, 

for i in the upper slab A, j in the down slab A, 

J 4s , for k, I on the surface, 
J 4fe , for i, j, k, I inside slab A, 
J^in, for i, j, k, I on the interface, 



(2) 



(3) 



Q s , for i on the surface, 

fit, for i inside ferroelectric slabs, (4) 
Clin, for i ontheinterface, 
where Jb is the nearest-neighbor coupling constant inside slab A. J s is the nearest-neighbor 
pseudo-spin interaction constant on the surface of slab A. J in is the nearest neighbor 
pseudo-spin interaction constant on the interface, is the distance between lattice site 
i and j (r^ > 1), a is introduced to describe the magnitude of the long-range interaction. 
Here, the long-range interaction is supposed to decay as 1/r 3 at large distance [22], and 
we take a = 3 in this paper. Considering that the environment of different site at the 
same layer is identical, we assume that the average value of the pseudo spin in the same 
layer of the ferroelectric slab has the same value. Within the frame work of the mean-field 
approximation and the decoupling approximation < s z -s z k sl >~< s z - >< s z k >< sf > , the 
average pseudo spin along z direction in ith layer can be expressed as following [16]: 



R. =< Sf >-- 



— , 1 , tanh - — — , 
2HA 2k B T' 



where 



< H Z >— ^ JijRj + X] JijklRjRkRl + 2/xi?, 



(5) 



(6) 



jkl 



H + (< h z >y, 



(7) 



where % runs over all the sites within the ferroelectric multilayer structure. In order 
to make the computation practicable, the long-range interaction is cut off at the 14th- 
neighbor in our calculations. It indicates that the value Ri is related to Rj for \rij\ < 4 
(according to our cut-off approximation). The error of the cut-off approximation can 
be neglected and the feasibility of the cut-off approximation has been discussed in our 
previous work [17]. When i runs over all the ferroelectric layers in the multilayer, the above 
Eq. (5) forms a set of simultaneous nonlinear equations from which Ri can be calculated 
numerically. Then, the polarization of the ith ferroelectric layer is 

nuRi, for i G the two — sided ferroelectric slabs, 
P, = I (8) 

0, for % G the non — ferroelectric slab, 
where n is the number of pseudo spins in a unit volume. The average polarization of the 
multilayer is obtained: 

1 N 

P(IV = ~y ^] Pjl (9) 
1 = 1 

where N is the total layers of the multilayer structure(including the non-ferroelectric lay- 
ers). 



The dielectric susceptibility of the ferroelectric multilayer is defined as 

dP(E,T) i 1 " dR i (E,T) i 

X(T) = — ^— I** = ^^2^—^—1^,. (10) 

i=i 

The deviation ^\ E=Q can be obtained by numerical differential calculation, then x(T) is 
obtained numerically. By changing the position and the thickness of the non-ferroelectric 
slab, we investigate the effects of the non-ferroelectric slab on the polarization and the 
susceptibility of the ferroelectric film. 



III. RESULTS AND DISCUSSIONS 



As reported [9], the increase of the lattice parameter c is caused by the mismatch of 
the in-plane lattice parameters between BaTiO^j ' SrTiO^. Here we take Jj n > J& in order 
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to ensure the enhancement of the polarization on the interface at a finite temperature. 
The coupling constant J s on the surface is weaker than the bulk in the general case. In 
this paper, we select Jj n = 2.0, J s = 0.5, and = 1.0. Different ratios of the four-spin 
and two-spin interaction coupling constant in the pure ferroelectric bulk material and 
ferroelectric film will result in many interesting phenomena [18, 19], such as the change of 
the order of the phase transition. If J s /Jb < 1 and Ji s / J&b < 1, the polarization and the 
phase transition temperature are reduced at the surface. If Ji n / -h > 1 and Jn n / -hb > 1, 
both of the polarization and the phase transition will increase at the interface [15]. We 
take the four-spin interaction coupling constant as J± s /J s = J^jJb = J 4m/ Jin — 6.0. The 
selected parameters will ensure that the order of phase transition of the corresponding 
bulk ferroelectric material is first-order, and the polarization of the interface between the 
ferroelectric layer and the non-ferroelectric layer will be enhanced. 

Firstly, we investigate the effect of the thickness of the film on the polarization and 
the phase transition temperature of the film with fixed thickness of the middle non- 
ferroelectric slab. Let L b = 2, and L a = L' a = 2,5,8. From Fig. 2(a), we can see that 
the average polarization and the phase transition temperature of the film will increase 
with the increase of L a (L' a ). And the surface polarization also follows this tendency (See 
Fig. 2(b)). If L a = L' a , the ferroelectric film with a middle non- ferroelectric slab can be 
considered as two symmetric coupled ferroelectric film without defects. So, the increase 
of the thickness of the ferroelectric film with a fixed thickness of the non-ferroelcetric slab 
is equivalent to the increase of the thickness the ferroelectric film without defects. It is 
the reason why the polarization and the phase transition temperature increase with the 
increase of the thickness of the film with fixed value of L b . 

From Fig. 3, we can see as the increase of the non-ferroelectric slab, the average 
polarization of the film will decrease, and susceptibility of the film will also decrease. 
There are two reasons for the decrease of the average polarization. The first reason is 
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the increase of the component of the non-ferroelectric slab. The second reason is that 
the corelation between the upside slab and the underside slab becomes weak when the 
thickness of the non-ferroelectric slab increases. The decrease of the susceptibility is 
resulted by the increase of the component of the non-ferroelectric slab that is irrelevant 
to the external electric field. 

In order to investigate the position effects of the non-ferroelectric slab, the thickness 
of the non- ferroelectric and the total thickness of the film are fixed as Lb = 1, L tota i = 11. 
The effects of the position of the non-ferroelectric slab on the polarization of the film 
and the upper surface are given in Fig. 4(a) and 4(b) respectively. It is found that the 
polarization of the upper surface and the film decrease with the decrease of the distance 
between the non-ferroelectric layer and the upper surface at low temperature, and increase 
with the decrease of the distance at the region of high temperature. The phase transition 
temperature goes up with the decrease of the distance. From Fig. 4, we can see a step-like 
polarization, which is also found in the pure ferroelectric thin film within the framework 
of the Landau Theory [20]. As the shift of the position of the non- ferroelectric slab to 
the upper surface, the upper surface will become disordered more easily, compared with 
the whole film, as the increase of the temperature. Supposing only the nearest neighbor 
interactions are considered, the surface phase transition temperature will be lower than 
the film. This step-like polarization is not corresponding to the ferroelectric-ferroelectric 
phase transition defined in ref.[21]. Here no phase transition is taken place at the tem- 
perature (T s ) where a step-like polarization occur. Due to the exist of the long-range 
interaction, the upper surface is still a ferroelectric state for T > T s . We can safely say 
that the occurrence of the step-like polarization is associated with a partial disorder of 
the pseudo-spins at the surface. If the long-range interaction is not taken into account, 
the surface will become complete disordered at T = T s . The more the non- ferroelectric 
slab near the surface, the larger polarization "step" will occur. 

The change of the susceptibility with the change of the position of the non-ferroelectric 
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slab is given in Fig. 4(c). When the position of the non- ferroelectric slab shifts to the 
upper surface, two peaks is observed in the susceptibility curve, and the two peaks cor- 
respond to the disorder of the pseudo spins in the upside slab and the underside slab of 
the film respectively. As the decrease of L a , the first peak will shift to the region of low 
temperature, and the second peak will shift to the high temperature region. When the 
non- ferroelectric slab sites on the middle position of the film(L a = L' a — 5), only one peak 
of the susceptibility is found. This is because the upside slab and the underside slab are 
symmetric when the non-ferroelectric slab sites on the center of the thin film, the two 
peaks will merged into one. 

In summary, we have investigated the polarization and the susceptibility of a ferroelec- 
tric film with first-order phase transition. The effects of the non-ferroelectric slab on the 
film can be concluded as follows: (1) The increase of the thickness of the non-ferroelectric 
slab will lower the polarization of the film and the susceptibility of the film. (2)When 
the position of the non-ferroelectric shifts to the upper surface, a step-like polarization is 
found in the surface, and the phase transition temperature of the ferroelectric film will 
be heighten, the number of the peaks of the susceptibility will change from one to two. 
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CAPTIONS OF FIGURES 

Fig. 1: 

The schematic graph of a ferroelectric multilayer structure with a non-ferroelectric slab. 
Fig. 2: 

The surface polarization and the average polarization of the ferroelectric multilayer versus 
temperature for different thickness of the ferroelectric slab L a (L' a = L a , L b — 2). 

Fig. 3: 

The surface polarization, the average polarization, and the susceptibility of the ferroelec- 
tric multilayer versus temperature for different thickness of the non-ferroelectric slab L& 
(L' a = L a = 5). 

Fig. 4: 

The surface polarization, the average polarization, and the susceptibility of the ferroelec- 
tric multilayer versus temperature for selected positions of the non-ferroelectric slab in 
the ferroelectric film. L b = 1, L tota i = 11. 
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